####################################################################################
###
###   Immigration, Public Housing and Support for the French National Front
###   Gloria Gennaro
###
###   Paper Table A16
###
####################################################################################


library(dplyr)

################################################################################
# Set Up
################################################################################

# Load Working Directories
source(paste0(wd_main, '/2_code/00_working_directories.R'))

# Load data and clean
source(paste0(wd_code, '/01_data_load_clean.R'))


################################################################################
# Auxiliary Functions
################################################################################

df = df_did[which(df_did$year=="1999"), ]


vars = c("fn_1995", "imm_share_90", "imm_share_99", "pop_99", "median_inc_01", 
               "less_primary90", "highschool90", "asso_pc99", "hlm_share_90", "res_owned_share_90",
               "empty_log_share_90", "congestion", "popgrowth")


mean_table = aggregate(. ~ imm_quant_99, 
                        data = df_did[, c("imm_quant_99", vars)], 
                        FUN = function(x) round(mean(x, na.rm = TRUE), 3))  %>% t()

print(mean_table)
